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[57] ABSTRACT 

A method for determining Richardson Number, Ri, or 
its reciprocal, RRi, for clear air turbulence prediction 
using measured potential temperature profiles and de- 
termining the vertical gradient of potential temperature, 
d0/ dz. Wind vector from the aircraft's instrumentation 
versus potential temperature dw/d# is determined and 
multiplied by dO / dz to obtain dW/dz. Richardson 
Number or its reciprocal is then determined from the 
relationship 

for use in detecting a trend toward a threshold value for 
the prupose of predicing clear air turbulence. Other 
equations for this basic relationship are disclosed to- 
gether with the description of a method for combining 
of other atmospheric observables using multiple regres- 
sion techniques to produce improved warnings of clear 
air'turbulence. 

7 Claims, 8 Drawing Sheets 
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MICROWAVE TEMPERATURE PROFILER FOR 
CLEAR AIR TURBULENCE PREDICTION 

ORIGIN OF THE INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract, and is sub- 
ject to the provisions of Public Law 96-517 (35 USC 
202) in which the Contractor has elected not to retain 
title. 


TECHNICAL FIELD 

The invention relates to a method employed aboard 
an aircraft for predicting clear air turbulence (CAT) 
from remote sensor measurements of air temperature 
vertical gradient and in situ measurements of horizontal 
wind vector from an inertial navigation system (INS) 
for inferring vertical wind shear, and more particularly 
to determining from these measurements atmospheric 
dynamic stability (Richardson Number), which is gen- 
erally recognized as having CAT encounter predictive 
capability. 

BACKGROUND ART 

CAT continues to plague commercial and military 
aircraft. CAT represents more than just an annoyance 
to passengers; it is a distraction to the flight crew, it 
consumes extra fuel (by reducing fuel efficiency and by 
requiring horizontal diversions around regions of sus- 
pected CAT), and it contributes to airframe metal fa- 
tigue leading to structural failure during flight and land- 
ings. 

A new type of unmanned vehicle is being designed 
for military reconnaissance. Some of these vehicles 
have large wingspans, fly slowly, and cannot tolerate 
moderate turbulence because of their vulnerability to 
unrecoverable aerodynamic upset. These reconnais- 
sance aircraft also need a CAT prediction capability in 
order to avoid such turbulence. 

In the field of CAT prediction, key factors that deter- 
mine Richardson Number (Ri) are vertical air tempera- 
ture gradient and vertical wind shear. Since low values 
of the Richardson Number are known to be necessary 
for production of CAT, it is important to have on board 
an aircraft the facilities for monitoring Ri, or its recipro- 
cal RRi, for predicting in real time when and where 
CAT will occur and how severe it will be. 

STATEMENT OF THE INVENTION 

In accordance with the present invention, Richard- 
son Number (Ri) is determined for CAT predictions. 

There are four steps to providing the information 
needed for determining Ri. The first step is to determine 
the vertical gradient of potential temperature, d0/dz. 
This may be done by measuring brightness temperature 
with a passive microwave radiometer operating at fre- 
quencies near 58 GHz, where oxygen molecules emit 
thermal radiation. The viewing direction of the radiom- 
eter is scanned through a selection of elevation angles. 
Measured brightness temperature versus elevation 
angle is converted to air temperature T versus altitude, 
z. Potential temperature 0 is then derived from air tem- 
perature and pressure at the aircraft's altitude. These are 
then converted to profiles of potential temperature, 0, 
from which it is possible to derive the vertical gradient 
of potential temperature, d0/dz. 

The second step is to determine the correlation of 
wind vector (from the aircraft’s INS) versus potential 


2 

temperature 0, and more specifically by determining the 
slope of both the EW and NS wind components. U and 
V, versus potential temperature, i.e., dU/d0 and dV/d0. 
The orthogonal sum of these is the “gradient of wind 
5 vector with potential temperature," dW/d0. 

The third step is to multiply dW/d0 by the “gradient 
of potential temperature with altitude" to obtain "verti- 
cal wind shear," 

10 



The fourth step is to calculate Richardson Number, Ri, 
^ which is the ratio of stabilizing forces to overturning 
forces derived from the above-mentioned vertical gra- 
dients: 


20 




It is more convenient to use the reciprocal RRi— I /Ri 
as the parameter for CAT prediction. When RRi < 1 for 
an atmospheric layer, it is stable and the growth of small 
vertical perturbations is suppressed. When RRi >4 for 
an atmospheric layer, small perturbations grow until 
wave breakdown. Typically RRi is 10“ 1 to 10“ 3 , and 
only rarely is as large as 1. The trend of RRi along a 
flight path allows for extrapolation, and when the ex- 
trapolated RRi trend extends toward the critical region 
of 1 to 4, a CAT warning is indicated. 

The severity of the turbulence can be determined by 
35 the magnitude of overturning energy per unit altitude 
increment, which is simply an evaluation of (dW/dz) 2 
calculated for the extrapolated situation of RRi = 4, 
w r hich can be calculated from d0/dz and dW/dz. CAT 
avoidance in response to a warning calls for an evasive 
4 Q maneuver based on the premise that most CAT is pro- 
duced by Kelvin-Helmholtz wave growth within inver- 
sion layers. Therefore, the most severe CAT is to be 
avoided by making the quickest flight level change 
leading to escape from the inversion layer. The temper- 
45 ature profile can be used to locate the altitude bound- 
aries of the inversion layer. 

The novel features that are considered characteristic 
of this invention are set forth with particularity in the 
appended claims. The invention will best be understood 
50 from the following description when read in connection 
with the accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates a stack of air masses in the strato- 
55 sphere overlying each other with shallow transition 
layers between them. 

FIGS. 2a through 2d depict time traces of an imagi- 
nary flight through a transition region that has high RRi 
values during the middle portion of the flight. 

60 FIGS. 3 a through 3/ show three 20-second sets of 
wind tracks for both EW and NS horizontal compo- 
nents, U and V, as a function of potential temperature 0. 

FIG. 4 is a plot of VWS and lapse rate (dT/dz versus 
time). 

65 FIG. 5 is a plot of CAT intensity versus time during 
a flight segment. 

FIG. 6 is a plot showing Richardson Number, Ri, and 
CAT intensity during the flight segment of FIG. 5. 
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FIG. 7 is a plot of the reciprocal of Richardson Num- 
ber, RRi, and CAT intensity versus time for the same 
flight segment of FIGS. 5 and 6. 

FIG. 8 is a functional block diagram of an airborne 
microwave temperature profiler (MTP) for a CAT 
warning system. 

FIGS. 9 a and 9b illustrate examples of a cockpit dis- 
play that could be generated using the concepts of this 
invention. 

DETAILED DESCRIPTION OF THE 
INVENTION 

The stratosphere consists of a stack of air masses 
overlying each other, with shallow transition layers 
between them, as shown in FIG. 1. Each air mass moves 
in a distinct horizontal direction governed by its his- 
tory. Consequently, vertical wind shear is concentrated 
in the transition layers. Temperature contrasts are also 
concentrated in the transition layers which are often 
“inversion layers." The most common generator of 
CAT is the transition layer separating the troposphere 
from the stratosphere, called the tropopause. A tropo- 
pause containing an inversion layer temperature feature 
is an especially CAT prone transition layer. 

The presence of the transition layers can be inferred 
from vertical profiles of Richardson Number, or prefer- 
ably the reciprocal of Richardson Number, RRi, i.e., 
the ratio between overturning forces to stabilizing 
forces. The transitions are marked by enhanced values 
of RRi. Since air masses can have a component of verti- 
cal motion, the interfaces are not always horizontal. 
Flight at a constant altitude may sometimes intersect the 
transitions layers, as shown in FIG. 1. 

CAT is generated when Ri<j. i.e.. when RRi ex- 
ceeds 4, and it persists until RRi becomes less than 1. It 
is still open for speculation how far RRi will be driven 
below I by the homogenizing action of CAT. In non- 
CAT regions, RRi typically has values of about 0.1 or 
less. Since air for which RRi >4 will produce CAT 
throughout neighboring regions meeting the criterion 
of RRi > 1. there is always the possibility of CAT when 
RRi> 1. Further, since CAT is felt at layers that border 
the altitude regime where it is generated (albeit at a 
lower severity level , it is important to keep in mind that 
CAT can even be encountered when RRi<l. If RRi 
varies gradually with altitude, then larger RRi values 
can be expected to be associated with greater probabili- 
ties of CAT, even when RRi< 1. Thus, the probability 
of CAT at given intensity levels can be expected to 
increase in going from “typical" RRi to “borderline" 
RRi, where typical ~ 0.1 and borderlines 1. 

It is also important to keep in mind that CAT reduces 
RRi, bringing it closer to “typical" values, with a time- 
scale that may be short compared to the life of the CAT. 
In other words, CAT may quickly “w'ash out" the high 
RRi values that generated it, then linger at ever- 
decreasing severity, during which the energy in the 
vortices cascade to smaller and smaller scales (eventu- 
ally becoming “thermalized"). Hence, if CAT produces 
RRi values that are <1, yet are still >0.1 (the normal 
background level), these regions should have a finite 
probability of being turbulent due to the difference in 
timescales for RRi destruction and CAT dissipation. A 
conclusion from these several arguments is that RRi 
values below 1 may not be as “innocent" as would be 
suggested by first-order Kelvin-Helmholtz wave theory 
and Richardson Number stability theory. 
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FIGS. 2a through 2d depict time traces of an imagi- 
nary flight through a transition region that has high RRi 
values during the middle portion of the flight. In FIG. 
2a, it can be imagined to be a horizontal transition layer 
5 that is close to instability for only the middle portion, or 
it can be imagined to be produced by flight through a 
transition layer that is slightly inclined to the horizontal 
and which is close to the threshold for becoming unsta- 
ble throughout a much larger horizontal extent than is 
10 showm. The following arguments w ill be valid for either 
imagined geometry. 

FIG. 2a show's RRi versus time, and the central part 
of this plot exhibits RRi above normal and slightly 
exceeding RRi = 4 at the center. In the early part of the 
15 flight, CAT has not yet been generated because the 
instability criterion RRi >4 has only been reached at 
location L5. 

FIG. 2b shows what might be observed on a later 
flight through the same region, after an outbreak of 
20 CAT has occurred at location L5 where RRi had ex- 
ceeded 4 for a sufficient time. FIG. 2c shows conditions 
during an even later flight pass, after several patches of 
CAT have occurred at locations L2-L6. At this time, 
the CAT from these outbreaks is still dissipating, which 
25 takes longer than the initial homogenizing of the transi- 
tion layer (and reduction of RRi). Thus, CAT may still 
be encountered where RRi has been reduced, but it is 
not encountered where RRi is still high (as well as 
where RRi has never been high). 

30 In FIG. 2d. all the high RRi regions centered on 
locations L:~L6have produced CAT, and the transition 
layer has “exhausted itself’ of CAT-producing vertical 
wind shear. All RRi values are <1, and perhaps they 
are much less than 1. However, residual CAT (at some 
35 intensity level) may still exist throughout the middle 
region of FIG. 2d. If the remnant RRi values distinguish 
themselves from the normal background levels, then the 
“intermediate" values for RRi might be associated with 
light CAT. 

40 Note that there is a negative correlation of CAT w'ith 
RRi within the region of initially high RRi. The dis- 
tance scale of the initially high RRi (the field of CAT 
patches) is typically 500 to 1000 km, whereas the dis- 
tance scale for the individual CAT patches is 20 to 50 
45 km. Therefore, with this model, the method of the pres- 
ent invention predicts a negative correlation of CAT 
with RRi for distance scales of less than 100 km, and 
positive correlations for distance scales exceeding 100 
500 km. 

50 Kelvin-Helmholtz (KH) wave amplitude is theorized 
to grow- unchecked and trigger CAT, provided RRi >4 
for a sufficient length of time. CAT sustains itself (i.e., 
continues to extract energy from the horizontal wind 
field and thermalize that energy) until RRi^l. When 
55 KH waves overturn, they homogenize the CAT layer 
by stochastically exchanging air parcels. This process 
leads to a reduction of the vertical wind shear to very 
low- values and also produces a fairly uniform potential 
temperature field (thereby completely destroying any 
60 pre-existing inversion layer). Thus, the method of the 
present invention will allow* for the possibility that 
CAT converts air from RRi ^ 4 to RRi < 1 . Perhaps RRi 
is reduced to values much lower than 1, such as 0.1. 

It should be kept in mind that RRi and CAT have 3-D 
65 distributions, and not 2-D as suggested by FIGS. 2 a-2d 
for simplicity of illustration. The sequence of CAT 
formation and spreading can be expected to have dis- 
tinct patterns in the vertical dimension. For example, 



5,117,689 


5 

CAT may star! at a middle altitude (where RRi is pre- 
sumably a maximum) and spread vertically toward both 
the upper and lower edges of the transition layer. An- 
other complication to keep in mind is that in the real 
world flight paths can course through these distribu- 5 
tions at odd angles and in other than straight and level 
flight. This complication will have to be dealt with in 
further refinements of CAT warning algorithms for 
ascents and descents. 

Vertical Wind Shear (VWS) Algorithm 10 

In the method of the present invention, determining 
VWS requires measurements of (1) the U and V wind 
vector components by an inertial navigation system 
(INS), where U and V are the EW and NS components 15 
of the horizontal wind vector, ( 2 ) air temperature and 
aircraft pressure altitude, and (3) air temperature lapse 
rate, defined to be dT/dz, derived by a microwave 
temperature profiler (MTP) shown In FIG. 8 . Air 
“brightness temperature" versus elevation angle is mea- 2 o 
sured by a vertically scanning radiometer and con- 
verted to air temperature T versus altitude z. Potential 
temperature 0 is then derived from air temperature T(z) 
and pressure altitude P(z) using the conventional defini- 
tion: 25 

0^T(K]<KXX)/P[mb3> o::8t5 . 
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independent variables. The proposed regression analy- 
sis produces solutions for the horizontal gradient of the 
wind components, dU/dt and dV/dt. the gradients of 
the wind components with 0 , dU/d 0 and dV/d0. and 
the RMS fit before and after the fitting analysis, RMS 0 
and RMS/. 

Three-point weighted averages of the partial deriva- 
tives of U and V with 0 are less affected by horizontal 
wind gradients, provided a proper weighting equation is 
used. The weighting can be performed on the U and V 
gradients separately using the following equation: 

Weight = !/{(! +(HG/5.6)2 + (RMS/35) 2 } 

where HG = dU/dt or dV/dt, and has dimensions [m/s 
per second]. The weighted-average dU/d0 and dV/d0 
can be converted directly to dU/dz and dV/dz by mul- 
tiplying by d0/dz. When this is done, the orthogonal 
sum of the resultant derivatives of U and V with z cor- 
responds to the desired VWS. 

Another equivalent algorithm is also possible, which 
preserves like-observables for as long as possible in the 
derivation (and corresponds with the derivation in the 
section above). The orthogonal sum of the two 0 - 
derivatives of U and V can be formed, called 
“VWS0"[m/s per OK]: 


Surfaces of equal 0 are called “isentropes" and since 0 is rwse={(dU/de) 2 +(dl'/d0) 2 }l ( 1 ) 

a property of the air that does not change with an air 

parcel's adiabatic altitude excursions, isentropes can 30 By multiplying “VWS0" by “d0/dz," the value 
usually be thought of as streamlines of air movement. VWS [m/s per km] for the vertical gradient of horizon- 

The crucial insights leading to the realization that tal wind in terms of altitude z (instead of 0) is obtained. 

VWS can be derived with existing INS sensors, and d0/dz [0K/km] is calculated from MTP-derived lapse 

w hich is used in the present invention, is summarized by rate dT/dz [K/km], and in situ (INS) measurements of 

the following two facts: * air pressure, P [mb], according to the relation: 

( 1 ) dU/d0 and dV/d0 can be derived using in situ mea- 
surements of U and V while 0 is also derived; this is de/dz=(dT/dz+]0) m (\ooo/P) 0 - 2 * b (2 

done by creating scatter plots of U and V versus 6 

during altitude excursions of either the aircraft or By combining these equations, we obtain a relation for 
isentrope surfaces, i.e., by real-time correlations of 40 VWS [m/s per km]: 
wind vector components U and V with potential 

temperature 0 for specified intervals, such as 30- rH5= vwse ( de/dz ) (3) 

seconds. The 30-second data for the U and V wind 

vector components are regressed against 0 and time Ri is normally stated in terms of the parameters just 
simultaneously as independent variables to minimize 45 described, namely: 
the confounding influence of horizontal wind gradi- 
ents on the solution for the vertical gradient. Ri-(g/0) [(d6/dz)/vws 2 ] ( 4 ) 

(2) These 0-derivatives can be converted to altitude 

derivatives by simply multiplying by d6/dz, which where g = 9.8 [m/s 2 ], the force of gravity. However, a 
can be calculated from knowledge of air pressure and 50 simple rearrangement of terms allows an alternative 
air temperature lapse rate dT/dz (measured by the expression to be .derived for Reciprocal Richardson 
MTP instrument). Number, RRi. That expression may be regarded as 

These two steps allow calculation of dU/dz and dV/dz, better because it contains quantities that are “closer" to 

which can be orthogonally added to produce the vector the observed quantities. The expression is: 

value of vertical wind shear, VWS, which is the crucial 55 

parameter needed for calculating Richardson Number. */?/- ( 0 / 9800 ) (dT/dz+ \0) { 1000 /p ) 0 286 

A more compact equation for Richardson Number (and {(dU/dB) +(dV/dO) } 

its reciprocal) is derived in the following section appen- 
dix, and it is this equation which is used in the analysis where 0 has units [KJ, dT/dz has units [K/km], and 
of the present invention. 60 dU/d0 and dV/d0 have units [m/s per 0K], The recip- 

It will be instructive to follow derivations of more rocal Richardson Number (»1/Ri) is better “behaved" 
familiar atmospheric properties. The next several para- I^ an Ri w'hen temporal averaging is performed. Conse- 

graphs are presented with this purpose in mind. quently, RRi is preferable to use for the derivation of 

For various reasons there may be an optimum dura- GAT warnings, 

tion for the time chunks of data used in performing least 65 FIGS. 3 a to 3/illustrate three 20 -second sets of wind 

squares regression analyses in which the U and V wind versus potential temperature tracks that “jump" indicat- 

components are chosen, in turn, to be the dependent ing a vertical w'ind shear, VWS. FIGS. 3a c and e are 

variable, and time and 0 are used simultaneously as for EW-component (U) data and FIGS. 3 b. d and /are 
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for NS-component (V) data. Note the pattern of lower- 
left to upper-right for the left column plots and the 
upper-left to lower-right slope for the right column 
plots. The EW-component U of VWS is positive, while 
the NS-component V of VWS is negative. The total 5 
VWS, being the orthogonal sum of U and V compo- 
nents. is approximately 30[m/s per km] during this per- 
iod. 

The left column of FIGS. 3 a, c and e shows that for 
three successive 20-second intervals the U component io 
increases with 0, with a slope of about -f 1.0 [m/s per 
0K], and the V component decreases with 0 with a 
slope of —2.0 [m/s per OK]. For this case, the gradient 
of the wind vector shears at the rate of 

15 

[m/s per BK\. 

FIG. 4 shows dT/dz averaging - 1.5 [K/km] at a time 
corresponding to the data of FIG. 3 (time equals 63.500 
seconds). Adopting 0 = 468 K and P 56 [mb], RRi is 20 
calculated 

RRi= 0.93 Udlvde) 2 ^(dV/d0) 2 ) 

where the 6 wind gradient is in units [m/s per 0K]. For 
a V WS0 gradient of 2.2 [m/s per 0K], calculated above, 
it is calculated that RRi = 2.0. or Ri = £. which is close 
to the critical threshold for Kelvin-Helmholtz instabil- 
ity considering that RRi typically has values <0.1 (i.e., 

Ri is typically > 10). 

30 

Example Illustration of RRi and CAT 

FIG. 5 portrays CAT intensity versus time during the 
same flight path depicted in FIG. 4 in which horizontal 
flight path position is given in terms of time, or ”UT 
seconds” with respect to an arbitrary reference time. 35 

FIG. 4 plots lapse rate dT/dz and VWS versus time 
FIG. 6 is a plot showing Richardson Number. Ri dips 
below 2 for the first time at 62.5 ksec (62,500 sec past 
the reference time) and 8 minutes later, at 63.0 ksec, 
light CAT is encountered. There is a pattern of low Ri 40 
from 62.5 to 69.0 ksec, where Ri is generally <5. Out- 
side this region Ri is typically zr 10. The solid line trace 
for Ri is actually the reciprocal of running averages of 
RRi. which accounts for the trace tending to follow the 
lower region of Ri excursions plotted in dotted lines 45 
instead of the upper region. CAT intensity is plotted 
using the scale on the right. 

Note that in the plot of FIG. 6, there is a correlation 
between CAT spikes and Richardson Number ap- 
proaching and dropping below 1. This is made much 50 
clearer by plotting the reciprocal of Richardson Num- 
ber, RRi, in FIG. 7. The reason for the better behavior 
of RRi is related to the different way fluctuations in the 
values of numerators and divisors affect their ratio. It is 
always better for a parameter to be defined with the 55 
most fluctuating part as the numerator. VWS exhibits 
greater percentage fluctuations in the real atmosphere 
than d0/dz, and it appears squared in the Ri equation. 
Since Ri has VWS in the denominator, it is preferable to 
plot its reciprocal, RRi, for display and use it in predict- 60 
ing CAT. 

The 2.5 minute averaged RRi trace reaches a maxi- 
mum of 1.9 at 63.4 ksec. The unaveraged data shows a 
maximum of 2.2 at this time, consistent with the data of 
FIG. 3 , which was for this time. Immediately after the 65 
peak in RRi, one of the main CAT events begins. There 
is a negative correlation of CAT with RRi within the 
early part of the “CAT containing region.” CAT 
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patches extend from 63.0 to 67.0 ksec (referring to FIG. 
5 ). The pattern of high RRi extends from 62.5 to 69.0 
ksec, a slightly larger region. 

It is possible to imagine that before the outbreak of 
CAT patches, the RRi trace rose from 0.1 to values > 1 
and perhaps as high as 4 throughout the “CAT contain- 
ing region.” As each CAT patch occurred, it quickly 
reduced RRi to low values and was left surrounded 
(horizontally) by high RRi regions. The edges of the 
high RRi region are assumed to be the last to generate 
CAT; i.e., that CAT patches occur first near the center 
of the high RRi region and then randomly occur at 
locations that move toward the outer edge of the high 
RRi region. At some stage in this process, there may 
exist an outer ring of high RRi, inside of which RRi is 
low but turbulent (with decaying CAT). This assump- 
tion w'ould provide an ideal geometry for predicting 
CAT encounters. 

The RRi criterion for distinguishing between CAT 
patch association and no association seems to be 0.2 (i.e., 
Ri = 5). Within the “CAT containing region,” which is 
defined as including the “ring” regions (i.e., from 62.5 
to 67.6 ksec), RRi exceeds the proposed RRi criterion 
859£ of the time. Outside the “CAT containing region” 
RRi exceeds the RRi criterion only 169c of the time. 

Measuring “vertical wind shear” in this manner and 
thereby inferring Ri or RRi, has not been previously 
known or suggested by others. It is a new tool that 
should, on theoretical grounds, have relevance to CAT 
prediction. There may be errors in estimating VWS, 
leading to VWS underestimates or possibly overesti- 
mates. It is nearly impossible to verify that the VWS 
values are correct, short of flying at least two aircraft in 
formation, one slightly above the other. 

If the RRi inferred from the measured VWS contin- 
ues to correlate with CAT experiences, the question of 
whether the new VWS rneasurements are accurate be- 
comes moot for the user interested in predicting and 
avoiding CAT. The motto “If it works, use it!” would 
then apply. 

Multiple Regression Combination of Observables 

Prediction of CAT likelihood should make use of 
more observables than just Ri or the reciprocal RRi and 
the trend thereof. From the model depicted in FIGS. 
2a-2d, which is supported by the data of FIG. 7, the 
correlation of CAT with RRi can be complicated. 
However, “multiple regression” analysis techniques can 
cope with this situation provided an archive of real 
measurements of observables and CAT intensity is 
available for analysis. 

An archive of observed parameter values, such as 
RRi or d0/dz, in conjunction with subsequently- 
measured near-future CAT intensity, can always be 
analyzed using multiple regression analysis techniques 
to yield a set of linear multiplying “retrieval coeffici- 
ents” which enable the conversion of any single set of 
parameter values to a near-future (predicted) CAT in- 
tensity. Past CAT itself can be considered an “observ- 
able” in this context. The trend of past CAT is another 
observable. The fluctuation level of air temperature in 
the recent past and the trend of this parameter are other 
“observables.” There are a potential 20 or so airborne 
measurable properties that merit consideration for this 
analysis. Clearly, some of these are more relevant than 
others, and the multiple regression analysis will assign 
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larger coefficients to these observables compared with 
the others. 

Ri or the reciprocal RRi and the trend thereof, are 
extremely likely to be one of the most important, most 
highly correlated observables for near-future CAT pre- 5 
diction. Until now it has not been possible to measure Ri 
or the reciprocal RRi. Analyses without Ri or the recip- 
rocal RRi show useful correlations with near-future 
CAT, and these other observables therefor contain 
information which should not be ignored in making 10 
CAT predictions. It is therefore very desirable that a 
multiple regression scheme be used to make use of all 
observables that can be measured besides Ri or its recip- 
rocal. Thus, the previous argument for how Ri or the 
reciprocal RRi can be used is meant only to illustrate 15 
the rationale for the concept that Ri or the reciprocal 
RRi have predictive value in forecasting CAT. There- 
fore, Ri or RRi should be included in an observed data 
base of atmospheric conditions on which a “multiple 
regression*' analysis technique is used to formulate for- 20 
casts of CAT. 

Implementation of Measurements 

A passive micro wave radiometer is used to measure 
brightness temperature while scanning the radiometer 25 
through selected elevation angles, such as within the 
range — 50 c to +60\ Measured brightness temperature 
at the selected elevation angles is converted to air tem- 
perature T as a function of altitude, z. This air tempera- 
ture data is then converted to profiles of potential tern- 30 
perature 0 versus altitude. The slope of 6( z), defined as 
d#/dz, is one of the required parameters needed in de- 
riving Ri or RRi. An airborne microwave temperature 
profiler (MTP) adaptable to a CAT warning system is 
shown in FIG. 8 and is described in U.S. Pat. No. 35 
4,346,595 and in a paper by Richard F. Denning. Steven 
L. Guidero. Gary S. Parks and Bruce L. Gary, “Instru- 
ment Description of the Airborne Microwave Temper- 
ature Profiler," J. Geophys. Res., 94, D14, 

16.757-16.765. Nov. 30, 1989. By reference to this paper, 40 
the instrument description contained therein is hereby 
made a part hereof. 

The MTP is comprised of two major assemblies, a 
sensor unit 10 and a data unit 11. shown in FIG. 8. A 
superheterodyne radiometer 10a is a major component 45 
of the sensor unit which is installed in an aircraft and 
used to measure profiles of air temperature versus alti- 
tude. It operates at frequencies near 58 GHz where 
oxygen molecules emit thermal radiation. Brightness 
temperature is measured at a selection of ten viewing 50 
elevation angles every 14 seconds in order to provide 
observed real-time quantities that contain information 
about the temperature of the air along the viewing di- 
rection, which is normally selected to be in a vertical 
plane in the direction of flight path but may also be in a 55 
vertical plane to one side. The observed quantities are 
converted to air temperature versus altitude z by the 
data unit 11 using pressure altitude of the aircraft as the 
origin or reference for converting viewing elevation 
angles to altitude. 60 

The altitude coverage by the ten viewing elevation 
angles from —50° to +60° is from about 2 km below to 
3 km above the aircraft altitude for flight altitudes of 
60,000 feet or 18.3 km, which were typical for the previ- 
ous usage of this particular MTP instrument. By chang- 65 
ing observing frequency slightly, the same 5 km region 
(from —2 to +3 km) can be sampled at lower altitudes. 

An MTP that observes at many frequencies in the 55 to 
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64 GHz region is able to provide 5 km altitude coverage 
for all cruise flight altitudes. 

Air temperature estimates are made for altitudes 
within this 5 km altitude range. Temperature profiles 
are obtained for small increments along the flight path, 
typically every 2 km for jet aircraft.. A super hetero- 
dyne radiometer 10c in the sensor unit is equipped with 
a horn antenna 106 and a reflector assembly 10c which 
may be adjusted in position by a stepper motor lOd to 
enable the reception pattern of the horn antenna to be 
viewing directions within the range —50° to +60°. 

The superheterodyne radiometer 10c is conventional. 
It employs a single mixer and two local oscillators (LO), 
one operating at 57.3 GHz and the other at 58.8 GHz, 
which correspond to a wavelength of about 6 mm. Four 
or more LO units can be incorporated to allow for 
operation at more altitudes. A microprocessor 11c in 
the data unit turns on one local oscillator at a time en- 
abling the measurement of atmospheric brightness tem- 
perature at many frequencies for each viewing direc- 
tion. 

The power entering the horn 106 is linearly propor- 
tional to the “brightness temperature" in the viewing 
direction. The conversion of radio frequency power 
level to brightness temperature is achieved using a mi- 
croprocessor 11c in the data unit 11, which applies the 
calibration offset of a target 10c to calibrate the output 
level of the radiometer for a known brightness tempera- 
ture. Both system gain and calibration target offset are 
determined periodically in flight for this conversion. 

The offset calibration is obtained by viewing an ambi- 
ent target 10c that is at a known physical temperature. 
The target used should have a high microwave emissiv- 
ity, which assures that its brightness temperature can be 
equated with its physical temperature. A thermistor 
temperature sensor is mounted in the target, and the 
target is allowed to change temperature as the environ- 
ment temperature changes. This arrangement reduces 
the difference between the temperature of the target 
and the air outside the aircraft, which in turn reduces 
the influence of any gain inaccuracies. 

The gain calibration is obtained by injecting a known 
amount of extra microwave energy in the waveguide 
immediately behind the horn antenna. This calibration 
signal is generated by a broadband noise diode 10/ 
which has been calibrated before installation in the 
aircraft using external hot and ambient calibration tar- 
gets (described below). The noise diode exhibits a stabil- 
ity that is adequate for establishing system gain. 

The microprocessor 11 a in the data unit 11 controls 
the stepper motor lGd for scanning the sensor unit. The 
microprocessor takes readings of the radiometer output 
and various MTP temperatures and records these data. 
The data unit is connected to the aircraft instrumenta- 
tion 116 which includes an inertial navigation system 
(INS) to provide in situ EW and NS wind vector (U and 
V) information and other aircraft sensor inputs for stabi- 
lizing the radiometer viewing direction in the presence 
of aircraft roll and pitch and to monitor aircraft altitude 
and air temperature. 

When the MTP is turned on from the cockpit, the 
data unit 11 automatically begins a series of 14-second 
observing cycles. A cycle consists of measuring the 
radiometer output from both channels for a sequence of 
viewing directions. There are 10 “sky" elevation angle 
views, and one view directed at the ambient calibration 
target. At the calibration target position, the calibration 
noise diode is also turned on to establish the radiometer 
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gain for both channels. Note that the calibration noise of 
the diode 10/ is coupled into the waveguide output of 
the horn antenna 106. The temperatures of the radiome- 
ter window and diode are transmitted to the data unit 11 
as analog temperatures. 

Radiation received through the waveguide coupling 
of the horn antenna through an isolator lOg into the 
heterodyne radiometer 10 j is mixed in a mixer lOg with 
the alternately selected local oscillator operated at 57.3 
and 58.8 GHz. The oscillators are coupled to the mixer 
lOg by isolators 106 and 10/. At the output of the mixer 
lOg. an IF amplifier and detector 10/ produce a voltage 
signal proportional to IF amplitude (i.e., atmospheric 
brightness temperature). A voltage-to-frequency con- 
verter 10A converts the voltage signal to a pulse train 
the frequency of which is proportional to signal volt- 
age. The digital pulse train is transmitted to a micro- 
processor input unit 11c which receives that digital 
pulse train and other signals from the aircraft instrumen- 
tation including pitch, roll, U, V, T and P, the latter 
being pressure altitude. These are individually con- 
verted to digital form and stored until accessed by the 
microprocessor 11a The analog temperatures are trans- 
mitted in analog form to the data unit 11 where an 
analog-to-digital converter 11 d converts the analog 
temperatures in the various channels into digital form 
and stores them until accessed by the microprocessor 
11a A stepper motor controller 11c is driven by the 
microprocessor lid to scan the radiometer in elevation 
angle through the requisite ten positions plus one posi- 
tion for the calibration target 10c. 

The microprocessor 11<7 is programmed to carry out 
all of the data processing and analysis and to display on 
a CRT display ll/either selected data, such as altitude 
temperature profile. RRi or vertical wind shear (VWS), 
or text messages summarizing a warning of CAT likeli- 
hood. severity and recommended avoidance maneuver. 

FIGS. 9 a and 9b illustrate examples of a cockpit dis- 
play that could be generated using the concepts of this 
invention. FIG. 9 a is an ATP (altitude temperature 
profile) display which could be selected by the pilot. It 
is more likely that a message display will normally be 
presented. An example of a message display is shown in 
FIG. 96. The message display includes the following 
information: “TR" means that a tropopause transition 
layer is nearby (in altitude); the “ + 0.5" states that the 
base altitude of this feature is 5 FL units above the 
aircraft, where one FL (flight level) unit is 100 feet; the 
up-arrow is a suggested avoidance maneuver; the verti- 
cal stack of numbers “2, 43, 4” signify that a minimal 
avoidance could be achieved by making a flight level 
change of 2 FL units upward, or 4 FL units downward, 
and that if no flight level change is made there is a 43% 
chance of moderate or greater CAT during the next 10 
minutes. 


12 

-continued 
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40 I claim: 

1. A method for determining Richardson Number, 
Ri, or its reciprocal, RRi, for clear air turbulence pre- 
diction aboard an aircraft in flight comprising the steps 
of 

45 determining vertical gradient of potential tempera- 
ture, 66/ dz, 

determining the correlation of wind vector from said 
aircraft's instrumentation versus potential tempera- 
ture, d\V/d0, 

50 multiplying the gradient of potential temperature 
with altitude, d0/dz, by said gradient of wind vec- 
tor with potential temperature, dW/d0, to obtain 
vertical wind shear, dW/dz, and 
calculating Richardson Number, Ri, from said verti- 

55 cal gradients from an equation 


APPENDIX 
Derivation of Ri or RRi Equation 
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60 

or its reciprocal, RRi, for use in detecting a trend 
toward a threshold value for the purpose of pre- 
dicting clear air turbulence. 

2. A method as defined by claim 1 wherein said verti- 
65 cal gradient of potential temperature, d0/dz, is deter- 
mined by measuring brightness temperature with a pas- 
sive microwave radiometer operating at a frequency 
where oxygen molecules emit thermal radiation and 
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scanning said radiometer through selected elevation 
angles above and below the flight line of said aircraft, 
and converting measured brightness temperature versus 
elevation angle to air temperature. T, versus altitude, z, 
and then converting said profiles of air temperature to 
profiles of potential temperature. 0 , and from said pro- 
files deriving said vertical gradient of potential tempera- 
ture dO/dz. 

3. A method as defined in claim 2 wherein determin- 
ing the correlation of wind vector from said aircraft's 
instrumentation versus potential temperature, dW /dO y 
is accomplished by determining the slope of both the 
east-west and north-south wind components, U and V, 
versus potential temperature 6 to obtain dU/d0 and 
dV/dO and computing the orthogonal sum 

[(du/di?) : -(dv/de) 2 ]i 

to obtain the gradient of wind vector with potential 
temperature, d\V/d0. 

4. A method for enhancing the quality of dW/d0 in 
claim 3 by reducing the influence of horizontal gradi- 
ents of wind on the d\V/d0 solution by performing a 
simultaneous regression of wind vector components U 
and V against both 6 and time, so that the time correla- 
tion has minimum influence on the 6 correlation. 

5. A method as defined in claim 1 wherein Richard- 
son Number, Ri. is the parameter of clear air turbulence 
prediction when the trend of Ri extends toward the 
critical region of 1 to J from the normal range of 5 to 
100 or the reciprocal Richardson Number, RRi, is de- 
termined as the parameter of clear air turbulence pre- 
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diction when the trend of RRi extends toward the criti- 
cal region of 1 to 4 from the normal range of 0.01 to 0.2. 

6. A method as defined in claim 5 wherein the sever- 
ity of the clear air turbulence predicted is determined by 

5 evaluation of (dW/dz) 2 calculated for the extrapolated 
situation of RRi = 4, which is to say that the value of 
(dW/dz) 2 when RRi =4 can be calculated from the 
values of dO/dz and dW/dz during the time that RRi is 
trending toward the critical region \ to 4. 

7. A method for making predictions of near-future 
clear air turbulence aboard an aircraft in flight by com- 
bining many atmospheric observables and weighting 
each of said observables with an assigned coefficient 

15 previously determined by multiple regression analysis 
of a data base of the intensity of clear air turbulence 
recorded in the past together with attendant atmo- 
spheric observables, in which Richardson Number, Ri, 
or its reciprocal, RRi, is determined for clear air turbu- 
20 lences recorded in said data base from information in 
said data base for inclusion as an observable in said 
multiple regression analysis to determine a correspond- 
ing coefficient assigned thereto, and combining each 
observable measured in flight multiplied by said as- 
25 signed coefficient determined by multiple regression 
analysis of said data base, including as an observable 
said Richardson Number, or its reciprocal, determined 
in accordance with claim 1 multiplied by said corre- 
sponding coefficient determined by multiple regression 
analysis of said data base. 

* * ♦ * * 
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